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ABSTRACT 

Using the DR12 public release of APOGEE data, we show that the thin and thick disks 
separate very well in the space defined by [a/Fe], [Fe/H] and [C/N]. Thick disk giants 
have both higher [C/N] and higher [a/Fe] than do thin disk stars with similar [Fe/H]. 
We deduce that the thick disk is composed of lower mass stars than the thin disk. 
Considering the fact that at a given metallicity there is a one-to-one relation between 
stellar mass and age, we are then able to infer the chronology of disk formation. Both 
the thick and the thin disks - defined by [a/Fe] - converge in their dependence on 
[C/N] and [C+N/Fe] at [Fe/H]«-0.7. We conclude that 1) the majority of thick disk 
stars formed earlier than did the thin disk stars 2) the formation histories of the thin 
and thick disks diverged early on, even when the [Fe/H] abundances are similar 3) 
that the star formation rate in the thin disk has been lower than in the thick disk, at 
all metallicities. Although these general conclusions remain robust, we also show that 
current stellar evolution models cannot reproduce the observed C/N ratios for thick 
disk stars. Unexpectedly, reduced or inhibited canonical extra-mixing is very common 
in field stars. While subject to abundance calibration zero point uncertainties, this 
implies a strong dependence of non-canonical extra-mixing along the red giant branch 
on the initial composition of the star and in particular on the a-elemental abundance. 
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1 INTRODUCTION 

It is now well-established from analysis of many disk galax¬ 
ies that Milky Way-like galaxy disks have a two-component 
structure. A comprehensive set of 21 review articles pro¬ 
viding an overview of the development and present status 


of this knowledge is available in Oswalt & Gilmore (2013), 


with a detailed history of the evidence from star count stud¬ 


ies in particular available in Yoshii (2013). Modern stud¬ 


ies of the Galactic thick disk began with Yoshii (1982) and 
Gilmore & Reid (1983) utilising photographic star counts. 


Later, quantitative high-quality stellar spectroscopic studies 
showed that the local thin and thick disk have also distinct 
(^-element abundance enhancements: 


cf 


Reddy et al. 


Fuhrmann (1998); 


Adibekyan et al. 


(2003); Bensby et al. (2005) 

^2011); Recio-Blanco ( 2014| . Suggestions from analysis of 
the low-resolution SEGUE [a/Fe] data that the thick disk 


was not chemically discrete from the thin disk (Bovy et al. 


2012), in disagreement with both other analyses of SEGUE 


data (Lee et al. 2011) have been shown inconsistent with bet¬ 


ter resolution data in several later studies (Guiglion et al. 
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2015). We show below the very clear distinction between 


thin and thick disks manifest in APOGEE data. The sev¬ 
eral current high-resolution, high-quality large sample stellar 


spectroscopic surveys underway - Gaia-ESO (Gilmore et al 


2012]), APOGEE pdltzman fe et al.||2015) , GALAH flDe 
Silva &; et al.|2015 ) - together with the imminent Gaia data, 


are producing large precise data sets, to extend our knowl¬ 
edge far beyond either small samples or low-resolution ex¬ 
ploratory surveys. Recent chemo-dynamic models interpret 
available data in the framework of Galactic formation, with 
several different scenarios for Galactic disk evolution pre¬ 


sented, as for example the study of Haywood et al. (2013). 


Rix & Bovy (2013) provide a helpful recent review. A funda¬ 


mental challenge in interpretation of the rich available data 
in terms of a Galactic evolutionary history is the difficulty of 
assigning ages (relative or absolute) to field stars. Determin¬ 
ing stellar ages from observables will always require stellar 
evolution models, while stellar properties such as mass and 
luminosity are currently best constrained with asteroseis- 
mology and parallax measurements. While we await Gaia 
for mapping large samples of stellar populations, samples 
available for accurate age determination are currently lim¬ 
ited to a maximum of a few thousand stars. 






















































































2 T. Masseron G. Gilmore 


Low-mass giant stars (M < 3M©) are excellent tracers 
to probe the Galaxy in multiple dimensions (space, kine¬ 
matics and chemistry) because they are very numerous and 
they are intrinsically very bright. Therefore they can be ob¬ 
served and studied in detail up to very large distances in the 
Galaxy (and beyond) with modern instrumentation, includ¬ 
ing in particular large spectroscopic surveys such as Gaia- 
ESO and the APOGEE survey. 

The APOGEE survey obtains high quality spectra of 
Galactic giant stars in the H band, with a reduction sys¬ 
tem providing homogenously-determined stellar parameters 


and elemental abundances. The recent public release (Holtz- 


man & et al. 2015) provides a homogenous set of stellar 


parameters and elemental abundances, with in particular 
of relevance for this study, carbon, nitrogen and a-element 
abundances for more than 200,000 stars. While the calibra¬ 
tion of abundance and parameter scales derived from au¬ 
tomated mass-production procedures is a major study in 
its own right, with much progress anticipated by compar¬ 
ing the various surveys underway with each other, and with 
Gaia and asteroseismology results over the next few years, 
already with appropriate care this large dataset provides a 
powerful opportunity to study Galactic star formation his¬ 
tories, and both stellar and Galactic evolution. 

Of particular interest among the elemental abundances 
which can be derived from H-band spectroscopy are carbon 
and nitrogen. The C and N elemental abundances in red 
giants are sensitive to stellar characteristics including mass 
and age. Carbon and nitrogen participate in many nucle- 
osynthetic reactions so that their abundances are altered 
by all stars during their evolution. During some stellar evo¬ 
lutionary phases these altered abundances are dredged up 
into the stellar atmosphere, providing direct evidence of the 
stars nucleosynthetic history, which probes stellar age and 
mass. One must therefore disentangle stellar evolution from 
Galactic evolution, but carbon and nitrogen provide the in¬ 
formation to constrain both star formation and galactic evo¬ 
lutionary histories. This analysis can be further extended 
by utilising the a-element abundance, which is unchanged 
during stellar evolution, and is a known indicator of star- 
formation histories, and distinguishes the thick and the thin 
Galactic disks. 

In the present paper, we take advantage of the mea¬ 
surements of C, N and the (^-elements in the large sample 
of giants in the APOGEE survey and consider their dis¬ 
tribution in the context of state-of-the art stellar evolution 
models. 


2 DATA SELECTION AND STELLAR 
EVOLUTIONARY STATE 
CATEGORIZATION 

We are interested in the subset of stars which have well- 
determined relative abundances for carbon, nitrogen, iron, 
and the alpha elements. In order to ensure optimal data 
quality, and given the very large APOGEE data set, we 
are able to sub-select from the full sample only those with 
signal-to-noise ratio >100, while retaining a statistically use¬ 


T e ff < 4000K. Therefore we also remove those cooler stars 
from our discussion: in practise such cool stars are in any 
case not relevant to this study, as the discussion below ex¬ 
plains. This leaves a sample of ~81000 stars having C and 
N measurements and which met these selection criteria. 

The resulting subsample is presented in Fig. ^ which 
illustrates the relationship between the stellar [C/N] abun¬ 
dance ratio, the stellar evolutionary state, via logg, with a 
further indication of [a/Fe] ratio, divided into two regimes, 
high ([a/Fe]> 0.15) and low ([a/Fe]< 0.15). This division is 
consistent with the definitions observed in other APOGEE 
studies ( |Nidever &; et ah |2014| |Anders 2014) and corre¬ 
sponds to the Galactic thick disk and the Galactic thin disk 
respectively (cf Fig. [4). The distribution function observed is 
manifestly inconsistent with a mono-abundance population 
in the [a/Fe] - [Fe/H] plane, as already observed by other 
independent studies (e.g. Fuhrmann|2011| Adibekyan et al. 


2011 Lee et ah||2011| Bensby et al.||2014 ). In the figure we 

can observe the different evolutionary status of the stars and 
the relationship with the C/N ratio. Since the pioneering 
work of Iben (1965), standard stellar evolution theory pre- 


ful sample size. In Holtzman & et al. (2015) the analysis 


team caution that they still have some challenges in de¬ 
termining accurate parameters for cooler stars, those with 


diets that after leaving the main sequence stars evolve up 
the sub-giant (subG) branch. It is expected up to this point 
that the surface abundances are unchanged since formation 
for normal single stars. Things change, however, as the star 
starts its evolution up the red giant branch (RGB). Model 
expectation is that the first dredge-up occurs, bringing some 
of the material synthesized during main sequence evolution 
up to the surface. In particular, this dredge-up enhances the 
nitrogen surface abundance at the expense of carbon. This 
expectation is supported by Fig. [T] where we observe an 
obvious systematic decrease of the [C/N] ratio during and 
after the first dredge-up phase. Moreover, it is well estab¬ 
lished that the amount of newly-synthesised material mixed 
to the surface is a function of stellar main-sequence mass. 
Hence, it is expected that the more massive is a star, the 
deeper is the dredge up, with more of the nuclear products 
brought to the surface. In other words, the more massive is 
a star the lower will be its post first dredge-up C/N ratio. 
Note also that there is a second-order dependence of the 
depth of the first dredge-up on stellar metallicity (see |Cha r- 
bonnel (1994) for a fuller discussion). After first dredge-up, 
as the star evolves along the RGB, more mixing occurs, de¬ 
creasing further the C/N ratio ( |Gratton et al. 2000; |Ma rteII| 
et al.|[2008 ). This observational evidence for non-canonical 
RGB mixing, that is a process that is not driven by convec¬ 
tion, is consistent with the thermohaline mixing mechanism 
for changes in elemental composition with stellar luminosity 
(Charbonnel &; Zahn 2007b). 

To take into account the several effects of stellar evo¬ 
lution on C/N ratios we divide the sample into its five 
main evolutionary stages (see Fig. [I]): the sub-giant branch 
(”subG”) for which we assume no or little change of the at¬ 
mospheric elemental composition since the time of star for¬ 
mation; the first dredge-up phase (”DUP”); the lower red 
giant branch (”low RGB”) - which shows the atmospheric 
elemental abundance composition as modified by the first 
dredge-up; the red clump luminosity region of the RGB, in¬ 
cluding the red clump itself (’’RGB + Clump”); and the 
upper RGB (’’upper RGB”) - where more CN processing 
from non-canonical processes may occur. 

Note also that for the following discussion we avoid the 
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first dredge-up range because it contains a mix of stars at 
different stages of their first dredge-up. We also avoid the 
region defined by 2.4 < log# <3.1 because there is a con¬ 
tamination by stars belonging to the red clump. Visual in¬ 
spection of the figure in the range of the clump ("RGB + 
Clump”) suggests that the C/N ratio range is larger than 
for the upper RGB. This would be surprising, since the C/N 
ratio is not expected to change after the RGB tip has been 
reached. This appearance is not real, but is an artefact due 
to 2 factors: i) the number of sample stars is much larger in 
the clump than in the upper RGB, which visually enhances 
the scatter on the plot; ii) the region contains an uneven mix 
of RGB stars, with relatively high C/N ratios, and clump 
stars, with lower C/N ratios. 

The discussion above demonstrates that [C/N] ratios in 
the red giants of the Galaxy observed by APOGEE are pre¬ 
dominantly determined by the star’s evolutionary state, and 
that we can use this ratio to infer the distribution of stellar 
properties, notably their relative masses. Although our dis¬ 
cussion primarily is based on a star to star comparison, it 
relies on two main requirements: i) The observed C/N ratios 
are internally on a consistent scale and do not suffer from 
any dominant bias correlated with stellar parameters ii) the 
C/N ratio of most stars at the time of their formation was 
close to Solar, [C/N]=0. We now demonstrate how and when 
those requirements can be justified and validated. 


2.1 Measurement random errors 


According to Holtzman & et al. (2015), random uncer¬ 


tainties on C and N abundances are respectively 0.06 and 
0.07dex. Systematic uncertainties will also be present - see 
below - but do not restrict the differential analysis we im¬ 
plement for the following reasons: 

• We restrict our analysis into specific evolutionary 
stages. This minimises a bias with stellar parameters by re¬ 
stricting to samples that have similar Teff and log g. Regard¬ 
ing metallicity, APOGEE data cover a large range. While 
our discussions remains insensitive to bias when comparing 
data at a given metallicity, we rely on the careful checks and 
calibrations from the APOGEE team and assume that the 
broad trends as a function of metallicity are real. 

• We select only good quality data (i.e. data with SNR 
> 100) to ensure that the measurement of line intensities is 
homogenous for all stars with similar stellar parameters. 

• We compute the star to star dispersion for 
[C/Fe] at a given metallicity and we obtained 0.08, 
0.07, and 0.05 for subgiants, lower RGB and upper 
RGB stars of the thin disk, while we obtained 0.10, 
0.09, and 0.06 respectively for the thick disk stars 
. Those values can be compared to the average of 
individual errors on C abundance of 0.05. We con¬ 
clude that the total scatter in [C/Fe] at any [Fe/H] 
and selected evolutionary stage is comparable to the 
quoted observational errors. As illustrated in Fig. [3j 
similar results are seen for [N/Fe] 

• It is recognized that errors in stellar atmosphere mod¬ 
elling (such as using ID vs 3D/hydrodynamical models) can 
be more pronounced at lower metallicity (e.g.Magic et al. 


(2013)). C and N are measured in APOGEE via molecu¬ 


lar features (respectively CO and CN). Therefore, using the 
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Figure 1 . (top panel) [C/N] vs logg for the selected sample of 
stars from the APOGEE survey with high signal-to-noise spec¬ 
tra, and with T e ff > 4000K. The stellar sample is colour-coded 
according to [a/Fe], with red points having [a/Fe] < 0.15 (pre¬ 
dominately thin disk) and black points having [a/Fe] > 0.15 
(predominately thick disk). Our empirical stellar evolutionary 
phases, correlated with stellar surface gravity, are indicated by 
vertical dashed lines. The stages are sub-giants (”subG”), dredge- 
up phase (”DUP”), lower red giant branch (’’lowRGB”), red gi¬ 
ant branch and red clump (”RGB+Clump”) and upper red giant 
branch (”upper RGB”). The lower panel repeats the stellar evo¬ 
lutionary indications, and in addition shows predicted theoretical 
[C/N] vs logg relations for three different masses and metallic- 
ities, as indicated, with the models being from |Lagarde et al.| 
|2012 ). The predicted locations of the red clump are shown by 
asterisks. The consistency between the data and the model trends 
is discussed in the text. 



fact that C and N abundances rely on lines forming in the 
same region of the atmosphere (the coolest part), taking 
their ratio should then naturally minimise any metallicity 
(and other stellar parameters) effects due to error in spec¬ 
tral modelling. 

To address in particular this last statement, we 
have completed some tests using the molecular 
features which are used for the measurement of C, 
N and O in the APOGEE data (namely OH, CO, 
and NH). For these, we compare the equivalent 
widths of each single line of those molecules for a 
given typical parameter of the APOGEE sample 
(T e ff s= 4500K,logg s= 2.5, and solar abundances) and 
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compare them with those obtained applying the 
main sources of astrophysical parameter random 
errors relevant to C and N measurements, i. e. 
T e ff + 250K, log g + 0.5, and[0/Fe] +0.1, with this latter 
source of uncertainty being especially relevant for 
molecular equilibrium concerns. In the upper panel 
of Fig. [2j we can observe that changes in stellar pa¬ 
rameters change the equivalent widths of molecular 
lines as expected. While we refer the reader to the 
work of Tsuji (1973) for a deeper understanding 
of the reasons for those changes, the net result is 
that both CN and CO molecules show very similar 
changes in their lines strength as a function of 
temperature and surface gravity changes. However, 
changes in surface gravity affect differently the OH 
and the CN molecules. Moreover, as illustrated in 
the third panel of the same figure, O abundance 
changes also affect OH and CN. In contrast, it very 
little affects the CO molecule strength, because 
we consider here C/O < 1 , so that the number of 
CO molecules is bounded by the number of free C 
atoms, and not of O atoms. To consistently evaluate 
the propagation of all those errors, we thus compute 
the changes in abundances in the bottom panel of 
Fig. [2j considering that O, C, and N abundance 
measurements are derived respectively from OH, 
CO and CN. It is remarkable to note in that figure 
that all the changes in parameters result in changes 
in abundances in the same direction and similar 
amplitude. Therefore, we conclude that systematic 
measurement errors on C and N abundances are 
expected to partially cancel out when considering 
their ratio. 


One may wonder what would be the net effect 
for other combinations of parameters changes. Actu¬ 
ally, we have demonstrated that the O abundance is 
directly or indirectly the main driver for differential 
changes between C- and N- bearing molecules, and 
thus to change in C/N ratios. Therefore, another set 
of parameters will lead to similar results. 

For this test, we have chosen the worst case scenario, 
meaning that the errors used in are significantly 
larger than those published in [Holtzman & et al.| 
(2015), and further that those random errors have 
been added up, while in a more realistic treatment 
random errors should be randomly distributed, thus 
reducing further the amplitude of any variations in 
abundances simulated here. 

Including random errors related to noise in a more 
realistic test will only affect the line to line disper¬ 
sion, and thus is not expected to change the average 
abundances. 


Based on all those arguments, the precision of the C/N ratio 
is at least as good as O.ldex at fixed [Fe/H] and evolution¬ 
ary state, which is sufficient precision to derive quantitative 
properties of the Galactic disks, which we do in the next 
section. 
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Figure 2. Upper panels: Differences in line equivalent widths of 
molecules present in the APOGEE spectra obtained when chang¬ 
ing stellar parameters from T e ff = 4500K,logg = 2.0, [O/Fe] = 
0.0 to the indicated values. Lower panels: Differences in C, N and 
O abundances measured when changing simultaneously all the 
stellar parameters. The overall change in parameters affects the 
abundances with the same sign and with comparable amplitude, 
supporting the idea that using those elements as ratios is robust. 


2.2 Measurement systematic errors 


The discussion above, which supports the published esti¬ 
mates for precision, does not hold for zero-point accuracy. 

Comparison of the [a/Fe] vs [Fe/H] distributions in 
APOGEE (Fig. [4| with those of Gaia-ESO (Mikolaitis & 
et al. ■ 20141 esp their Fig 3) suggests small systematic scale 
differences at the O.ldex level in element ratio and/or Fe- 
abundance zero point. These do not affect the discussion 
here. 

Checking the calibration of C and N abundances is 
harder, since, as the discussion above explains, these change 
during a star’s evolution. However, there is available a rela¬ 
tively small sample of sub-giants. These stars have not yet 
experienced dredge-up, so their surface abundances reflect 
the abundance of the inter-stellar medium from which they 
formed. Fig.[3]shows the ratios of [C/Fe], [N/Fe] and [C/N] 
vs [Fe/H] for the sub-giant sample. One’s expectation for 
the appearance of these figures is a function of one’s as¬ 
trophysical preconceptions and preferred chemical evolution 
model. The [C/Fe] vs [Fe/H] relation decreases slowly and 
smoothly with increasing metallicity, passing close to So- 
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Figure 3. Testing data scatter and zero points with pre-dredge- 
up subgiants. The three panels show [C/Fe], [N/Fe] and [C/N] as a 
function of metallicity for the APOGEE subgiant stars. There is a 
continuous decrease of C/N ratio with metallicity, consistent with 
plausible expectations from Galactic chemical evolution. Both low 
[a/Fe] (thin disk) and high [a/Fe] (thick disk) sub-giant stars have 
consistent values in the region of [Fe/H] where both populations 
are present. The definition of thick and thin disks is described in 
the text. We also note that [N/Fe], and consequently [C/N] is not 
equal to 0 at solar metallicity, suggesting a zero-point offset in N 
abundances. 


lar at Solar [Fe/H]. This is astrophysically plausible and al¬ 
ready observed by Nissen et al. (2014), suggesting the carbon 
abundance calibration is not unreasonable. The [N/Fe] vs 
[Fe/H] relation is roughly constant below Solar, with a sys¬ 
tematic increase above solar metallicity. This evolution looks 
plausible as one can expect the nucleosynthetic products of 
low mass stellar evolution, in particular N production, to be¬ 
come increasingly significant compared to supernovae prod¬ 
ucts, which have low N production, with increasing time. We 
do not of course know that the super-solar stars are indeed 
relatively young, so this remains a consistency statement, 
not a proof. As a striking example, Fig. [3] clearly shows that 
the [N/Fe] value at solar metallicity does not match the So¬ 
lar (0) value for the bulk of the stars, as already observed by 


Israelian et al. (2004), highlighting the preliminary calibra¬ 


tions for this element in the current APOGEE data. There 
is a zero-point offset in the [C/N] ratio of approximately 
0.2dex, apparently dominated by the offset in [N/Fe]. Al¬ 
though this does not affect the conclusions of this study, 
which are based on purely differential comparisons, we ap¬ 
ply an constant offset correction in the N abundance of +0.2 
for ah the APOGEE sample for the rest of the paper. 

The [C/N] ratio varies smoothly between low metallicity 
sub-giants and high metallicity sub-giants. This range, if 
a calibration artefact, could affect this study. However, we 
note that the [C/N] ratio at given metallicity is uncorrelated 
with [Fe/H] (lowest panel of Fig. [I]), at least according to 
models, which suggests that the trend is not a calibration 
error, and in any case will not affect our conclusions below, 
which depend entirely on the different behaviour of [C/N] 
vs [Fe/H] as a function of [a/Fe]. More importantly, the 
change in [C/N] after the stars have undergone their first 
dredge up is several times larger than the scatter in the sub¬ 
giant range, and is an order of magnitude larger than any 
systematic dependence of [C/N] on [a/Fe]. In summary, for 
the purposes of this study, we may be confident that the 
apparent variation of the [C/N] ratio in the sample is not 
a calibration artefact, but is dominated by internal nuclear 
processes and stellar evolution. 

For convenience, in the following discussions and figures 
we recalibrate the APOGEE data by applying an offset of 
+0.2 in N abundances to match the solar abundance. This 
offset has no effect on any of the discussion below unless 
noted, and we do not consider our adopted value of +0.2dex 
a precise or recommended number. 


3 THE THICK AND THIN DISKS ARE 
DISCRETE IN [a/FE] VS [FE/H] 

It has been established by many studies that the a-elements 
provide a discrimination between the thick and the thin 
disk, with the thick disk having systematically higher [a/Fe] 
abundance rations than does the thin disk at given [Fe/H] 
(Mikolaitis & et al. (2014)). Although questioned by |Bovy| 
et al. ~( 2012) in their analysis of low-resolution SDSS-SEGUE 


data, the present high-resolution SDSS-APOGEE data show 
a very clear distinction, with strikingly discrete distribu¬ 
tions. This is shown in Fig. [4] The distribution is into two 
broad maxima, the upper of which we identify as the thick 
disk, the lower as the thin disk. We emphasise that, for 
present purposes, the contours chosen are indicative. We 
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Figure 4. (top panel) The distribution of [a/Fe] as a function 
of metallicity for APOGEE upper RGB stars. The distribution 
into two broad maxima, with a clear minimum in number density 
between them is very evident. We identify the upper distribution 
within the density contours as thick disk, and the lower distribu¬ 
tion within the contours as thin disk. In the bottom panel data for 
the upper RGB stars are color-coded in two bins: [C/N]< —0.3, 
and [C/N]> —0.3. This discrimination corresponds to a combina¬ 
tion of intrinsic stellar properties, especially main-sequence mass, 
and stellar evolution, dredge-up and other mixing events, and not 
to Galactic spatial distribution or kinematics. The [C/N] selec¬ 
tion clearly correlates tightly with [a/Fe] confirming that the two 
identified stellar populations do form discrete populations. 
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make no attempt here to fine-tune the selection, with in 
particular the extension of the high-a group all the way up 
to [Fe/H]=0 being suggestive rather than robust. Rather 
the whole point is that the distinction between the domi¬ 
nant two populations is so clear. This bimodal distinction is 
not an arbitrary density-labelling, as we show in the lower 
panel, which shows the high correlation between the [a/Fe]- 
based population definition and the [C/N] ratio. Rather the 
stars identified by the density contours form two distribu¬ 
tions which are clearly different in their stellar evolutionary 
histories, and in their present-day age distributions. We con¬ 
sider that evidence further in the remainder of this paper. 

An additional advantage of using the [C/N] ratio is that 
this new population indicator relies on molecular features, 
which can be observed in relatively low resolution spectra, 
unlike the case with a elements. Thus C and N provide an 
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Figure 5. [a/Fe] as a function of metallicity for lower RGB stars. 
The C/N ratios are color-coded as indicated on the scale bar. We 
can observe a C/N gradient within the thin disk stars. This is 
consistent with the dependence of the first dredge-up on mass 
and metallicity, implying, unsurprisingly, that more metal-rich 
thin disk stars are systematically younger than more metal-poor 
stars. See the text for a more detailed discussion. 
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excellent and robust spectroscopic complement to popula¬ 
tion allocation of thin disk and thick disk red giant stars. 

We do not here attempt to analyse population num¬ 
bers, or selection effects. Clearly at the metal-poor end 
the APOGEE selection bias is significant. In the region 
where both populations are evident, essentially for stars with 
—0.6 ^ [Fe/H] ^ —0.2, there is no obvious selection effect 
which can generate a tight correlation between two group of 
elements with very different properties. 

The distribution function observed is manifestly incon¬ 
sistent with a mono-abundance population in the [a/Fe] - 
[Fe/H] plane. We do note also that there does seem to be a 
real excess over random scatter of stars near solar metallicity 
with high [a/Fe] values. We discuss these more below. 


4 CARBON AND NITROGEN ALONG THE 
RGB: DETERMINING AGES FOR THE 
THICK AND THIN DISKS 

The a-elements retain a memory of the material from which 
a star has formed. Carbon and nitrogen, on the other hand, 
convey age-dependant information on the evolutionary his¬ 
tory of the individual star. We now exploit those distinct 
properties of stars to discuss simultaneously the history of 
the thin and thick disks, and the imprints of stellar evolu¬ 
tion. 


4.1 The lower RGB stars: Star formation histories 
and relative ages of the thin and thick disk 

4-1.1 Observations 

We consider first those stars classified in Fig. [l] as in the 
lower RGB phase of stellar evolution, after first dredge-up. 

In Fig. [5] which presents the [a/Fe] vs [Fe/H] relation 
for the lower RGB stellar sample, we have colour-coded the 
[C/N] values. The gradient in colours can then be trans¬ 
lated into a gradient in [C/N] which is clearly correlated 
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with metallicity in the thin disk stars. As noted above, and 
as illustrated in the lower panel of Fig. [T] the value of an ob¬ 
served [C/N] ratio can be due only to dependencies on two 
properties of the specific star: a metallicity dependence and 
a mass dependence. Moreover, a more quantitative examina¬ 
tion of figure 5 of |Cha rbonnel (1994) tells us that for stellar 
masses below ~ 3M 0 d 0 t it is stellar mass which dominates 
the changing [C/N] ratio, while for higher stellar masses it 
is the metallicity. 

Fig. 0 presents the distribution of [C/N] as a function 
of metallicity for lower RGB stars, with the [a/Fe] defini¬ 
tion selection colour-coded. When comparing thin and thick 
disks stars, the degeneracy between the mass and metal¬ 
licity effects on the C/N ratio can be broken where the two 
disks have an overlapping metallicity range. Therefore, when 
examining Fig. [6] one can understand that the thick disk 
contains in general lower mass stars than the thin disks. 
Considering the fact that at a given metallicity there is a 
one-to-one relation between stellar mass and age, we can 
conclude that C/N ratios represents a clean proxy for time 
and that the thick disk population formed earlier than the 
thin disk population. 


4-1.2 Simulation of the [C/N] ratio in lower RGB stars 


To verify and quantify to some extent the evolution of [C/N] 
ratios with time, we employ here stellar evolutionary models. 
It is not trivial to reliably disentangle the effect of mass and 
metallicity on the C/N ratio to determine ages. Hence, to 
evaluate more quantitatively the relative ages of the thin 
and thick disk, we used the stellar model grid of Lagarde 
et al. (2012) to simulate the evolution of [C/N] ratio with 
time (Fig. [7). We first linearly interpolate the Lagarde et al. 


(2012) grid in time and mass with steps of respectively dt = 


10 years and dM = 0.02M©. Then, we derive the density 
of stars p as a function of time and C/N ratio using the 
following relation: 

rt+6t rCN+6[C/N] rlog gmax 

p{t , [C/N]) = / M~ a dMdCNdt 

Jt JCN 1 /log gmin 


where M is the initial star mass, t is the time, C and N 
are the carbon and nitrogen surface abundances at time £, 
and log# is the surface gravity as provided by the stellar 
evolution models. We also assume a Salpeter IMF for all 
the Galactic stars, thus a = 2.35, which is an acceptable 
approximation over the very small range in masses relevant 
here. We applied the same selection definitions based on 
log# as in Fig. [T] thus for low RGB stars log gmin = 2.7 and 
log gmax — 3.1. We chose St and S[C/N] to be respectively 
10 7 years and 0.01 dex as a compromise between numeri¬ 
cal limitations and simulation resolution. This resolution as 
well as linear oversampling of the initial grid, may lead to 
some uncertainties in the actual star densities, but in the 
following discussion only the relative evolution of [C/N] as 
a function of time matters. We also consider the impact of 
the initial [C/N] value on the post dredge up value. While 
the [C/N] value observed in subgiant stars at solar metal¬ 
licity is compatible with the value in the models, the [C/N] 
value in the subgiants at [Fe/H]=-0.54 is slightly enhanced, 
by 0.2 dex (see Fig. [ 3 ] but considering the zero point cor¬ 
rection in N) due to enhanced carbon at lower metallicities. 




time (Gyr) time (Gyr) 

Figure 8. Same as Fig[7]but with different logg cuts. In general 
the trends are not affected by changes in our selection criteria, 
which only affects the scatter. 


Hence, we add this value to our final stellar density value p 
when running the low metallicity models. This should be a 
reasonable approximation to mimic the dilution of the first 
dredge up products in slightly [C/N]-enhanced material. We 
consider the consistency between the model abundances and 
the observed values further below. 

As expected, there is definitely a one-to-one relation be¬ 
tween the value of [C/N] and age. What is also striking with 
this simulation is the tightness of the relation, making the 
C/N ratio potentially a very robust time indicator, provided 
the models are a good match to data. The slope of the re¬ 
lation in Fig. [T] is such as to allow age determinations with 
formal precision of order 1-2 Gyr. Hence we confirm that 
thick disk stars are older than thin disk stars - according to 
the models. 


4.1.3 Simulation uncertainties 

However, as illustrated in Fig. [l] there are discrepancies be¬ 
tween the models and the observations. Obvious examples 
evident by visual inspection include the predicted vs the ob¬ 
served location of the Red Clump, and the inconsistency be¬ 
tween the predicted low values for [C/N] at low log g and the 
observed high values. One can expect our simulation results 
to be affected by such disagreements, in particular regarding 
our selection criteria in logg. We do not know here whether 
the differences are a problem in the APOGEE calibration 
for surface gravities or whether the models are not accu¬ 
rate enough, and we further recall that we have offset the 
APOGEE [C/N] zero point to match Solar, and thus to im¬ 
prove agreement with the models. We can however test the 
model sensitivity, and so have run the same simulation with 
different logg cuts. We demonstrate in Fig. [8] that varying 
the logg selection does affect relative stellar densities, but 
that overall the trends remain identical, as do our conclu¬ 
sions. 
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Figure 6. C/N ratios as a function of metallicity for lower RGB stars (dots), colour coded based on the [a/Fe] population discrimination 
as thin disk stars (red) and thick disk stars (black). Density contours are also displayed to guide the eye. We recall that N abundances 
are zero-point corrected by +0.2 from the original data. The distributions of [C/N] ratios are different between the two populations, 
showing that they have different age and star formation histories. 
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Figure 7. Simulation of the evolution of [C/N] in an RGB star population with metallicity [Fe/H]=0 (left panel) and [Fe/H] =-0.54 
(right panel) as a function of time. The masses listed indicate the main sequence mass of the stars dominating in density at a given 
time and metallicity. The models adopt the stellar grid from[Lagarde et al. (j201^|. The red and the grey areas represent the values as 
observed in Fig. [6] in respectively the thin and thick disk. 


4-1-4 Simulation results for the lower RGB stars and 
discussion 


From Fig.[7]and using the results from our albeit preliminary 
stellar evolutionary simulations, we can now interpret the 
various trends observed for both disk populations. As seen 
in the simulations, the [C/N] value shows an overall decrease 
when the metallicity increases, mostly driven by the decrease 
of the initial [C/N]. This explains why the general trends for 
the thin and thick disk stars in Fig.[7]are negative. The fairly 
small range in [C/N] for the thick disk stars is significant. In 
fact, our simulation indicates that low metallicity thick disk 
stars have approximately the same age as do the most metal- 
rich thick disk stars. Indeed, this is a further proof that 
the thick disk had a fast star formation and self-enrichment 
history, consistent with the long-standing deduction from 
the enhanced [a/Fe] values. However, we observe for the 
thin disk stars a steeper trend in [C/N] as a function of 
[Fe/H] than for the thick disk. Interpreting the [C/N] values 


for the thin disk in terms of our simulations, this clearly 
indicates that there is an age difference of ~6Gyr between 
the metal-poor thin disk stars and the solar abundance thin 
disk stars. This implies that the star formation history has 
been much slower in the thin disk than in the thick disk. 


Interestingly, the thin and thick disk stellar trends seem 
to converge in their value of [C/N] at [Fe/H]«-0.7. This 
would imply that both populations have similar ages at this 
metallicity, which corresponds to an age of 10 to 12 Gyr, con¬ 


firming the suspicions of Haywood et al. (2013) and Bensby 


et al. (2014). Although the absolute time values are quite 


different, and thus shed light on uncertainties affecting mod¬ 
els and observations, our conclusions are qualitatively com¬ 
parable with the work of Haywood et al. (2013). Without 
detailed modelling of the sample selection function, spatial 
distribution, radial gradients, and so on, we cannot in this 
study deduce more specific age information. 

Although our analysis allows us to draw robust, albeit 
qualitative, conclusions when discussing the relative age and 
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Figure 9. C/N ratio as a function of metallicity for upper RGB 
stars, with thick and thin disk stars classified from Fig. ^]and with 
N abundances corrected by +0.2. While thin disk stars show a low 
and almost flat level of [C/N], the thick disk stars have a high 
[C/N] ratio. 


evolution of the thin and thick disk, we note that the ab¬ 
solute age derived from Fig. [7] for the thin and thick disk 
respectively, agrees satisfactorily with the standard expec¬ 
tations (and in particular the solar age). Note that we make 
important assumptions regarding the cross-comparison of 
abundance ratios in the thin and thick disk stars i) the abso¬ 
lute value of [C/N] in the APOGEE data analysis pipeline is 
not biased by the a abundance; ii) the first dredge-up occur¬ 
ring in the two disk stellar populations is identical, meaning 
that it is insensitive to the special properties of thick disk 
stars, particularly their high [a/Fe] abundances. Although 
any a-abundance dependence on stellar dredge-up does not 
appear as clearly in the current data as does the evidence 
for non-canonical extra mixing we discuss in the following 
section, this question requires further study. 


4.2 The upper RGB stars: constraints on 
non-canonical extra mixing? 

4-2.1 Observations 


We now consider the upper RGB stars. The [C/N] vs [Fe/H] 
distribution for upper RGB stars is shown in Fig [4] while the 
corresponding distribution of [C/N] with [Fe/H] is shown in 
Fig. [9] It is interesting to compare these distributions with 
those seen on the lower RGB, and with the model expecta¬ 
tions shown in Fig. [l] 

The differences between Fig(6] and Fig[9] are a signa¬ 
ture of stellar evolution. As discussed above, after the first 
dredge-up occurs, as a star evolves further up the RGB, 
its C/N ratio is further depleted by non-canonical extra 
mixing. The favoured physical mechanism responsible for 
this decrease is thermohaline mixing. Although this physi¬ 
cal mechanism is a well-known physical process (also known 


as salt finger effect and is observed in the oceans), Charbon- 


nel & Zahn (2007b) introduced it in their models and found 


that it is expected to set in after the first dredge-up when the 
star reaches the RGB luminosity bump, when the hydrogen¬ 
burning shell encounters the chemical discontinuity created 


inside the star by the convective envelope at its maximum 
extent during the first dredge-up. This supplementary mix¬ 
ing allows more material from H burning to be mixed and 
brought to the surface. More support for the idea that the 
driving mechanism is related to the RGB bump has been 
illustrated by the observation of Li abundance change along 
the RGB branch in a globular cluster ([Lin d et al.| |2009). 
Hence, because of this mixing, it appears natural that the 
C/N ratios are globally lower in Fig. [ 9 ] than in Fig[6] How¬ 
ever a closer inspection reveals that thin disk stars have a 
constant C/N ratio while thick disk stars show higher values 
and a trend. 


4-2.2 Simulation of the [C/N] ratio in upper RGB stars 


As for the first dredge-up it is expected that there is a 
dependence of the thermohaline mixing on the mass and 
metallicity of the star. To disentangle this degeneracy, we 
simulate the expected C/N using the same modelling ap¬ 
proach as we implemented earlier for the lower RGB stars 
(Fig. 10). The C/N ratio from the (corrected) observations 
broadly matches the model expectation for the thin disk 
stars, although the distribution of values tends to values of 
[C/N] more positive than [C/N] =—0.4, which is not con¬ 
sistent with the range predicted by the models. For thick 
disk stars, where the observed distribution function has 
—0.3 < [C/N] < —0.1, none of the models reproduces the 


observations. 

Since [C/N] evolution on the upper RGB is very sen¬ 
sitive to evolutionary state, one may be concerned by our 
selection criteria and the disagreement between the logg 
provided by the models and the observed values. This could 
be more crucial as our definition of upper RGB stars also 
encompasses several distinct stage of evolution (upper RGB 
branch, possibly some contamination by clump stars and 
even some early AGB stars). Once again to test the robust¬ 
ness of our conclusions, we run the simulation with different 
choices in logg (Fig. ED- It appears that, as for the lower 
RGB models, while log# selection affects the star densities 
in the diagram, the trends remain unchanged. 


4.2.3 Simulation results for the upper RGB stars and 
discussion 

In detail, when comparing the lower and upper RGB thick 
disk stars in Fig[6]and Fig [ 9 ] any decrease in C/N ratio along 
the RGB appears relatively weak, testifying that the non- 
canonical extra mixing has been significantly reduced. This 
is the first time that such an effect of reduced non-canonical 
extra mixing on the RGB branch has been observed in a 
large population. In the literature the only studied mecha¬ 
nism that could inhibit non-canonical extra mixing is a mag¬ 
netic field, but this is expected to affect only decendents of 
Ap stars (Charbonnel & Zahn 2007a), which represent only 
5% of the stellar population. Because the thin and thick 
disk star samples studied here have similar properties ex¬ 
cept a element composition, this suggests that this may be 
the important factor affecting non-canonical extra mixing. 
An extra hint towards this hypothesis is provided by Fig.[l2] 
where a C/N gradient proportional to the a enhancement 
is observed in the thin disk stars. On the other hand, at 
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Figure 10. Simulation of evolution of [C/N] as a function of time using stellar grid from Lagarde et al. (2012). The red and the grey 
area represent the values as observed in Fig. [ 9 ] in respectively the thin and thick disk. While the values of C/N ratios in the thin disk 
approximately match the expectation of stellar evolution models, it appears very difficult to reproduce the C/N ratios in the thick disk. 
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Figure 11. Same as Fig |10| but with different logg cuts. In general 
the trends are not affected by changes in our selection criteria, 
which only affects the scatter. 


Figure 12. [a/Fe] as a function of metallicity for upper RGB 
stars. The C/N ratios are color-coded. We can observe a C/N 
gradient proportional to a enhancement for the thin disks stars. 


lower metallicities and even higher a- enhancement, it has 
been well-established since the work of S pite et al.| (2005]) 
that canonical extra mixing is very efficient ([C/N] ~ —1.0) 
in a enhanced halo stars. If thermohaline mixing is in¬ 


deed the physical mechanism (Charbonnel &; Zahn 
2007b[ Lagarde et al.|2012 ), more studies may be re¬ 


quired to be able to evaluate its efficiency between 
very low metallicity halo stars and relatively more 
metal-rich thick disk stars, as well as the impact of 
various composition environments such as the one 
between thin and thick disks. 


5 PUTTING THE SAMPLE TOGETHER 
5.1 C+N vs O 

While a star’s C and N abundances are affected by the star’s 
evolution, the sum abundance, C+N, should be conserved 
regardless of the evolutionary status for a low-mass star, and 
is representative of its formation conditions. Using this prop¬ 
erty, we will in this section be able to constrain the formation 


of the Galactic population by now inferring the properties of 
the progenitors of the C+N. We summarise the distributions 
of various elements and element ratios for the APOGEE gi¬ 
ant sample in Fig. |13| We note without further discussion 
that the abundance scale of oxygen, through [O/Fe], is low 
compared to other recent studies, such as |Bertran de Lis| 


et al. (2015). Clearly our conclusions here are sensitive to 


further studies of calibrations. 

One can interpret the various trends in Fig.^Jin the fol¬ 
lowing way: [a/Fe] represents the contribution ratio of core 
collapse supernovae versus type la supernovae; [C+N/Fe] 
represents the enrichment of low and intermediate mass star 
AGB contribution versus massive star contribution. The life¬ 
times of low mass AGB stars to release their nucleosynthetic 
products is larger than that for SNIa, which are larger than 
for SN type II. 

However, Fig. [T3] shows that thick disk stars have 
formed with higher initial abundances of C+N than did thin 
disk stars, with this high abundance being indicative of the 
early stages of Galaxy evolution. One can also see that, while 
the C+N enhancement in the thick disk stars decreases only 












































C and N in the Galactic disks 11 


0.3 


0.2 

17 

0.1 


0 


-0.5 0 0.5 


[Fe/H] 



-0.5 0 0.5 

[Fe/H] 

Figure 13. Element ratios as indicated for all classified thick disk 
(high-a) and thin disk stars in the APOGEE RGB sample. Recall 
that we have offset the [N/H] zero point to match solar values. 



slowly with increasing [Fe/H], it decreases much faster with 
[Fe/H] in the thin disk. This is explained if the thick disk 
formed on a shorter time scale than did the thin disk, in 
the same range of [Fe/H], so that that SNIa had released 
more Fe in the thin disk. Although low-mass AGB stars are 
believed to be the major sources for C and N in the discs, 
we can still observe a decrease in [C+N/Fe], which indi¬ 
cates a larger contribution of Fe from SNIa. However, the 
slope of [C+N/Fe] vs [Fe/H] is lower than for [O/Fe], which 
means that the C+N/Fe decrease is partially compensated 


by the contribution of low and intermediate mass AGB stars. 


Bensby et al. (2005) measured [Ba/Eu] and found an early 


s-process contribution, which confirms that AGB stars con¬ 
tribute to the C+N enrichment of the thin disk early on. 
Moreover, we observe a ’’knee” in C+N at [Fe/H]^0 in the 
thin disk. If this additional increase in C+N is related to 
AGB stars, then this may mean that there is a significant 
change in the production of carbon for AGB stars above 
solar metallicity. 

The abundance of [C+N] is similar for both thin and 
thick disk stars at metallicity ~ —0.7. We have shown ear¬ 
lier another correspondence in [C/N] at a similar metallicity, 
indicating that thin and thick disk had similar ages at this 
point. Considering also the observations of |Bensby et al.| 
(2014), suggesting a lower limit for thin disk metallicity at 
—0.7, and the discovery of some kinematically selected thin 


disk alpha-rich stars at [Fe/H]« —0.7 by Haywood et al. 


(2013), we seem to seeing that, for stars now in the Solar 
neighbourhood, thick disk enrichment continued indepen¬ 
dently of the start of formation of the thin disk between 10 
and 12Gyr ago. 


5.2 The rest of the sample 

For this study we have empirically defined two dominant 
stellar populations, thick and thin disks, using as definition 
population density. We excluded from consideration, until 
now, sources outside the dominant groups. These apparent 
outliers will be a mix of errors, but also may represent both 
expected outliers, such as binary mass transfer stars, and po¬ 
tentially interesting sub-populations. These stars are identi¬ 
fied in the element ratio diagrams in Fig |14[ 


5.2.1 The Intermediate ol- enhancement population 


Fig. [14] shows the element ratio distributions of the stars 
which lie outside the dominant population groups. While the 
density of these stars is low, the majority do form a smooth 
extension of the trends seen in the main groups. Further 
study of measurement errors and scatter are required before 
a detailed discussion is feasible. Among topics of future in¬ 
terest are to determine if the data are able to make any use¬ 
ful statement as to whether the intermediate—a population 
at solar metallicity is a younger extension of a thick-disk like 
population, or is an older metal-poor thin disk (migrated?) 
population ( Adibekyan et al. (2011); Bensby et al.| ( [20T4| ), 
or something else. 


5.2.2 Super-solar population 

Potentially interesting sub-structure is evident at super so¬ 
lar metallicity in Fig. [9j suggesting that many stars with 
[Fe/H=+0.25 have [C/N]=—0.45, an unexpectedly high 
value. While this super solar population contains stars with 
intermediate C/N ratios, it does not show any particular a 
enhancement compared to the thin disk population (Fig. [9]). 
Note also that this population does not appear as clearly 
in lower RGB stars, as indicated in Fig. [6] This may sug¬ 
gest that we detect the effect of metallicity on non-canonical 
extra mixing on the upper RGB, rather than a true dis¬ 
tinct population of super-metal-rich old stars. In the latter 
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Figure 14. In the top panel we define the sample of stars which 
lie outside the two main groups defining thick disk and thin disk. 
The [C/N] (middle) and [C+N] (lower) distributions of these stars 
are presented in the other panels (including a zero point offset in 
N abundances). 
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6 CONCLUSIONS 


Defining the Galactic thick disk is an interesting challenge, 
with several combinations of kinematics and chemistry in 
use. We show here that there is a clear bimodality in the 
distribution function of stars in the [a/ Fe] vs [Fe/H] space. 
Further, we show that this bimodality is very highly corre¬ 
lated with stellar [C/N] ratio in red giant stars. The [C/N] 
ratio of an RGB star is determined by its evolutionary state, 
and in particular is a function of stellar mass, and hence 
age. Our observed correlation of properties allows us to de¬ 
termine the relative age distribution functions of the thick 
and the thin disk, showing that the thick disk red giants are 
indeed systematically older than the thin disk red giants at 
the same value of [Fe/H]. Thus the thick disk is older than 
the thin disk. The thick disk also has higher values of [a/Fe] 
ratios, implying also that it formed more rapidly than did 
the metal-poor stars of the thin disk. 

C/N ratio is an excellent and alternative probe for con¬ 
straining the relative masses of giant stars and hence their 
relative ages without the use of isochrones. As a quantita¬ 
tive comparison, determining an age of RGB branch stars 
from isochrone fitting with a precision of lGyr requires an 
accuracy in T e ff better than 50K while it requires an accu¬ 
racy in [C/N] of 0.05. It is very tempting to assign masses 
and ages to the whole sample using the C/N ratios. How¬ 
ever this exercise would first require us to prove that the 
absolute measurement of the C/N ratios in the RGB of the 
present sample is better that the actual excellent precision 
- and we emphasise that we see an offset of approximately 
0.2dex in the zero point of [N/Fe] for solar abundance stars 
in the current APOGEE dataset. 

We also stress that although there are several indepen¬ 
dent studies to derive stellar population ages (e.g. Haywood 
et al. ( 2013| )), all attempts rely on stellar evolution models. 


Unfortunately, despite the huge effort currently devoted, we 
show that the models still suffer from some caveats, perhaps 
even for long-studied effects such as the first dredge-up. In 
the particular case of red giant branch stars, we show that 
more detailed models of the first dredge-up and of canonical 
extra mixing is required, specifically their sensitivity to the 
initial C, N and a —element composition. 
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case, these stars may be related to the recent finding of Ko- 


rdopatis &; et al. (2014), revealing a local population likely 


to have radially migrated from inner regions of the Galaxy. 
This group may also be related to the a—enhanced popula¬ 
tion as seen in figure 9 of Nidever &; et al. (2014). However, 
more theoretical study of the impact of a—abundances and 
metallicity on the dredge up and on non-canonical mixing 
processes is needed before drawing any firm conclusions. 
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